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5h ' 1. Introduction 

Polarimetry has become a fundamental part of solar physics and is rapidly gaining impor- 
tance in many other fields of astrophysics. 1 One of the most challenging problems that one 
faces when doing astronomical polarimetry is the determination of the instrumental polar- 
ization introduced by the telescope in the observed beam. There are two approaches that 
have been used to deal with this problem. The first one is to place calibration optics at the 
telescope entrance, which allows one to feed light in a known state of polarization into the 
system. The measured Stokes parameters can then be used to determine the telescope re- 
sponse with a very high degree of accuracy. Two solar telescopes currently use this strategy: 
The Dunn Solar Telescope (DST), at the Sacramento Peak observatory (NM, USA) operated 
by the National Solar Observatories, and the German Vacuum Tower Telescope (VTT) at the 
Observatorio del Teide of the Instituto de Astrofisica de Canarias (Spain). The Swedish Solar 
Tower at the Observatorio del Roque de Los Muchachos (also of the Instituto de Astrofisica 
de Canarias) will soon have this capability, as well. 

The second approach to determine the telescope polarimetric response is to observe rela- 
tively well-known astronomical sources of polarization. This is sometimes done in night-time 
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astronomy by observing "standard stars". In solar physics, this may be accomplished with 
observations of a sunspot. The procedure 2 considers Milne-Eddington fits (MEF) to sunspot 
umbra spectra, from which the elements of the telescope Mueller matrix are adjusted to 
minimize the residuals of the MEF to the observed data. 

It is generally desirable to have a purely instrumental telescope calibration, so that no a 
priori assumptions on the observed targets are needed. Unfortunately, the construction of 
calibration optics for a large aperture (larger than ~1 m) is impractical. This had not been 
an issue in the past because existing solar telescopes have apertures that do not exceed 1 m 
in diameter, whereas larger night-time telescopes did not require absolute polarimetric accu- 
racies higher than what can be routinely achieved by means of standard stars. This situation, 
however, is rapidly changing due to the ever increasing accuracy demands of astronomical 
polarimetry and the development of new large-aperture solar telescopes, like the German 
telescope GREGOR 3 or the Advanced Technology Solar Telescope 45 (ATST). 

A previous paper 6 introduced the "beam-expansion" method (BEM) for the instrumental 
calibration of a large telescope. This work presents a different procedure, the sub-aperture 
method, which has some interesting advantages over the BEM (at least for the particular 
case of the ATST design). We shall hereafter consider the problem of calibrating the ATST, 
but the proposed procedure may be easily applicable to other telescopes as well. 

The sub-aperture method is based on the following basic idea. One starts by effectively 
turning the ATST into a small-aperture telescope. This is done by inserting a mask into 
the beam at any pupil image. The telescope sub-aperture is calibrated as usual by means 
of calibration polarization optics above the primary mirror. An astronomical polarization 
source (e.g., a sunspot) is observed through the sub-aperture and then again with the full 
aperture. The two observations are only a few seconds apart in time. The sub-aperture 
calibration is used to remove instrumental polarization in the first observation, which allows 
one to determine the actual Stokes vectors emerging from the astronomical source. Once 
these input vectors are known, the full aperture observation can be used to calibrate the 
telescope. 

2. Requirements 

The sub-aperture method proposed in this work requires the following elements. Note that, 
for the case of the ATST, requirement 1 is probably the most demanding one. 

1. A mechanical mount is needed to slide calibration optics in and out of the beam before 
the primary mirror. The calibration optics may have a small diameter (see next point) 
and consist of an achromatic linear polarizer and a retarder. The polarizer and retarder 
can rotate and be inserted and removed independently of one another. It is important 
to realize that the calibration optics does not need to be placed over the center of 
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the aperture. It could be at any location over the primary mirror, e.g. near the outer 
edge of the telescope aperture. This consideration probably simplifies the design of 
the mechanical mount, since it does not need to extend all the way into the center of 
the aperture. Moreover, there is an additional advantage for off-axis telescopes like the 
ATST. By having the sub-aperture closer to the axis of the system, rather than the 
aperture center, the instrumental polarization of the sub-aperture is minimized. This 
is of relevance because, as discussed in §4, the level of accuracy achieved is limited by 
the accuracy of the sub-aperture calibration. 

2. The diameter of the calibration optics (do) may be as small as desired, provided only 
that it fills an aperture large enough to observe the source of polarization. Typically 
this would be a sunspot. However, there may be times of low solar activity in which 
no sunspots are available. In such cases plage regions may be employed as sources. A 
suitable value for do is probably between 20 cm and 40 cm (see §4). 

3. A blocking mask that can be inserted in the beam independently of the calibration 
optics is required. The location of the mask can be any position in the system where 
an image of the telescope pupil is formed. It may be at the entrance (e.g., the dome 
shutter) or, perhaps more convenient, at any location downstream in a collimated beam 
(which would make it a small element). The purpose of the mask is to block the rays 
in the light beam that did not pass through the calibration optics. In other words, 
when the mask is in the beam the telescope aperture is reduced to do- The alignment 
between the mask and the calibration optics does not need to be very accurate. We 
simply require that only light that passed through the calibration optics can get through 
the mask. 

3. Calibration procedure 

The procedure proposed to calibrate the optical train is as follows: 

1. Slide calibration optics into the beam 

2. Slide mask into the beam. Only polarized light passes through. 

3. Calibration Stokes vectors are measured for various configurations of the calibration 
optics, in the same way as it is presently done with other solar telescopes. If (xt) 
represents the Mueller matrix of the small- aperture telescope+instrument system, we 
have that: 

S°J = (xt)St, . (1) 
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4. From (known) inputs S 1 ^ and (measured) outputs S°^, determine (xt). In general the 
measured Mueller matrix, that we denote by (xt)', is affected of small errors. 

(xt)' ~ (xt) . (2) 

5. Remove calibration optics, but leave mask in. 

6. Observe solar source of polarization (e.g., a sunspot) through mini-aperture of diameter 

Ssolar,! = ( X ^)^oZar,l ■ (3) 

7. Remove mask and repeat same observation using full aperture. Only a few seconds pass 
between the observations with and without the mask. If (XT) is the Mueller matrix 
of the full-aperture telescope+instrument system, we have: 

S°olar,2 = 0^T)S™ lar 2 . (4) 

8. Calibrate small-aperture solar observations S™{ arl to obtain S^ arl : 
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^solar,l — solar, 1 ~ ( X ^) $solar,l ■ (5) 

9. Use known S™ larl and measured S™{ ar2 to determine (XT). The following assumption 
is implicit: 

qin ^ qin / r>\ 

°solar,l — °solar,2 ■ \ v ) 

There are two important assumptions in the procedure as outlined above. The first one 
is given by Eq (2) and requires that the calibration of the sub-aperture must be done with 
sufficient accuracy. Possible errors in this step will affect the full aperture calibration through 
Eq (5). 

The second assumption is that of Eq (6). In order for it to be accurate, it is important 
that the two observations S so i m . t i and S so i art 2 be taken as close in time as possible. It is also 
important that the field of view remain the same, especially when doing slit spectroscopy. 
Possible discrepancies in the telescope pointing between the two observations need to be 
much smaller than the spatial structure of the polarization source. This would probably not 
be an issue when an adaptive optics system is employed to compensate for seeing-induced 
wave-front distortion. However, it is not yet clear that the adaptive optics system of the 
ATST will be able to operate while in calibration mode. 

The validity of these two assumptions, and the degree to which they limit the accuracy of 
the telescope calibration, is explored by means of numerical simulations in §4. 



4 



4. Numerical simulations 

4- A. The hound and hares tests 

This section describes some numerical experiments aimed at estimating the accuracy of the 
sub-aperture method and its limitations. The simulations consist of several "hound and 
hares" type of tests, in which one seeks to determine a telescope matrix which is known 
beforehand. I have used calibrated data from the ASP instrument to simulate a set of solar 
Stokes vectors {S l ™ larl and S™ lar2 )- The simulated solar data is observed through a small 
sub-aperture system (characterized by its Mueller matrix (xt)) and through the full system 
(with Mueller matrix (XT)). 

Obviously, the actual Mueller matrix of the ATST is still unknown. Using that of an 
existing telescope is probably not very helpful, since their Mueller matrices are likely to be 
very different. The approach taken in this work is to perform montecarlo simulations with a 
large number of randomly chosen Mueller matrices that represent the telescope. According 
to reference, 6 off-diagonal elements of up to ~5% are to be expected due to the off-axis 
geometry of the ATST. Therefore, it seems appropriate to consider matrices of the form: 

xt = 1 + 0.05 x a x 
XT = 1 + 0.05 x <7 2 , (7) 

where 1 is the 4x4 identity matrix and a\ is a particular realization of a 4x4 matrix whose 
elements are random real numbers obeying a normal distribution of width 1. 

Two different ASP maps are used in the simulations. The first one (Map A) 7 contains a 
sunspot and some plage regions. The second (Map B) 8 is a quiet Sun map and is used here 
to explore the feasibility of the sub-aperture method with quiet Sun signals. Fig 1 shows the 
observations and the three slit positions that have been used as solar polarization sources 
(sunspot, plage and network). 

Let us discuss the numerical experiments in some detail. We start by taking some spectra 
from one of the slit position marked in Fig 1. These will be our solar input Stokes vectors 
from the first observation \S™ lar X (A, y), where y denotes spatial position along the slit and 
A represents the spectral variable]. An observation through the telescope sub-aperture is 
simulated as: 

S:l r , 1 (\,y) = (xt)ST olar , 1 (\,y). (8) 

The variation of (xt) with wavelength is slow and can be neglected over the spectral range 
recorded by the detector (typically a few A). 

A second set of Stokes spectra S™ lar2 (\ y) is obtained from the same map, but at a slightly 
different slit position. Let us denote the offset between the two observations by n x . These 
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Fig. 1. Left: Continuum intensity maps. Right: Magnetograms (degree of 
circular polarization). Upper panels: Map A. Lower panels: Map B. Scales on 
horizontal and vertical axes are in arc-seconds. The vertical lines mark the 
three slit positions used in the simulations as representative of sunspot, plage 
and network, respectively. 
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vectors are observed using the full aperture: 

5rL, 2 (A,2/) = (XT)5™ ar>2 (A, 2/ ). (9) 

Random noise is added to both S™f arl and S™f ar2 . The noise has a normal distribution 
and has a larger amplitude for the case of S°^ ar x as a result of the smaller observing aperture 
(do). I have simulated an exposure time of ~4 seconds in both cases, which results in a noise 
amplitude of 5 x 10~ 4 /°fo (with do given in meters) for S°^ arl and 1.25 x 10~ 4 for S°^l ar2 . 
Obviously, it would be possible to use longer exposure times, thus reducing the noise in the 
observations. 

At this point we have the simulated solar observations and the telescope matrices used to 
generate them. Next we need to apply the procedure described in §3 to the simulated data 
and compare the retrieved matrix (XT)' to the actual (XT) generated earlier. 

In a real situation, one starts by using the calibration optics to obtain the Mueller matrix 
(or at least an approximation) of the sub-aperture system (xt). The matrix (xt)' resulting 
from the measurements should be a good approximation to the actual (xt). For instance, 
the ASP calibration yields an accuracy 9 better than 5 x 10 -4 . 

The tests presented in this section consider that: 

(xt)' = (xt) + Vt a , (10) 

where rjt is a parameter that accounts for the accuracy of the sub-aperture calibration (e.g., 
for ASP rjt < 5 x 10~ 4 , as stated above) and a is again a matrix of random elements. 
Several experiments with different values of rjt have been carried out in order to explore the 
propagation of errors from the calibration of the sub-aperture to the full aperture system 
(see below). 

Finally, one needs to perform steps 8 and 9 of the procedure detailed in §3. A least-squares 
algorithm is employed to determine the (XT) matrix by solving the following system: 

S™lar,2 = (X-T)S™ lar 2 . (11) 

Instead of using all the possible Stokes vectors observed at every (A, y) point along the slit, 
it is numerically convenient to select a suitable less redundant subset. The strength of the 

-* in 

polarization signal of the S' solar ± is used as a selection criterion. A total of 280 vectors, 
consisting of the 40 points with the largest positive values of Q, -Q, U, -U, V, -V plus the 40 
points with lowest total polarization, are selected. Although the number of vectors has been 
chosen rather arbitrarily, the procedure is not very sensitive to it. This value was adopted 
after verifying that, in selected simulations, doubling it resulted in no accuracy improvement. 

The least-squares fitting results in a (XT)' matrix that can be straightforwardly compared 
to (XT). In this manner it is possible to study the propagation of errors and the overall 
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accuracy of the sub-aperture method. The entire process has been repeated 100 times for 
each one of the simulations, with different realizations of the random matrices that appear 
in Eqs (7). The results presented below are the mean-square errors from the montecarlo 
simulations. 

4-B. Results 

Let us start the discussion in this section with a simple reference case. Consider for instance 
a diameter of d = 0.30 m for the calibration optics and a calibration accuracy r) t = 5 x 10~ 4 . 
The solar source is a sunspot (slit position 1 in Fig 1). Let us assume for the moment 
that the telescope pointing has not changed between the two observations (rj x = 0). The 
propagation of errors and noise through the calibration procedure (as explained above) leads 
to the following error matrix E = (XT)' — (XT): 

/ 6.9 x 10" 4 1.3 x 10~ 3 2.4 x 10" 2 7.1 x 10" 4 \ 



E 



6.0 x 10~ 4 1.4 x 10~ 3 2.3 x 10~ 2 7.2 x 10~ 4 
6.2 x 10~ 4 1.2 x 10~ 3 2.5 x 10~ 2 7.3 x 10~ 4 
^ 6.8 x 10" 4 1.3 x 10~ 3 2.3 x 10" 2 7.1 x 10" 4 j 



(12) 



Notice that the second and third columns, which correspond to cross-talk from Q and U, 
are larger than the other two. This is due to the fact that the linear polarization signals are 
considerably weaker than I and V, even in the penumbra where the linear polarization is 
strongest. Fortunately, the ATST polarimetric accuracy requirement (the amount of cross- 
talk must be smaller than 5 x 10 -4 ) does not necessarily imply that the matrix elements of 
E have to be all smaller than 5 x 10 -4 . Typically, V is ~ 10 _1 / c (where I c is the continuum 
intensity) and Q,U are ~ 10 _2 J C . This means that, in order to meet this requirement, the 
elements in the first column of E must be smaller than 5 x 10 -4 , those in the second and 
third smaller than 5 x 10 -2 , and those in the fourth smaller than 5 x 10~ 3 (see reference 10 
for a more detailed discussion). These values are used as derived requirements for the ATST 
polarimetric accuracy in the discussion below. 

Fig 2 shows how the calibration error varies with the sub-aperture diameter do for r\t = 
5 x 10~ 4 and rj x = 0. The horizontal dotted line marks the ATST accuracy requirements for 
each column of (XT). Cross-talk terms from I are almost at the required level for d > 20 cm 
(upper left panel). A slightly lower r\t would be needed to fully meet the ATST requirement 
(see below). The situation is better for the other Stokes parameters, but only when using a 
sunspot as the polarization source. Plage and quiet Sun fields can be used to determine the 
first and fourth columns of (XT) at ATST levels. However, the second and third columns 
(cross-talk from Q and U) exhibit very large inaccuracies. The reason for this is that the 
observations used here have been taken at disk center and show very little linear polarization. 
In order to fully determine the Q and U terms of the telescope matrix from plage or quiet 
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Sun regions, one would also need near-limb observations which exhibit linear polarization 
signals. 

It is possible to achieve the required accuracy by calibrating the sub-aperture to the level of 
rj t ~ 3 x 10~ 3 , as shown in Fig 3. Notice in that figure that the telescope calibration accuracy 
is strongly related to that of the sub-aperture. Also in Fig 3 we can see the uncertainties 
introduced by pointing errors between the observation of S so i arj i and S so i ar ^- The method is 
not compromised by errors as large as 0.75 arc-seconds and is therefore robust enough to be 
employed even in the absence of adaptive optics. 

5. Advantages and Disadvantages 

The sub-aperture method has the following advantages over the BEM and MEF: 

• Unlike MEF, no a priori assumptions are made on the observed targets (physical model, 
profile shape or symmetry, etc). 

• Simpler design modifications than BEM. 

• Unlike MEF, can be employed at any observable wavelength. 

• Unlike BEM, uses the exact same optical setup as the actual observations. 

• Unlike BEM, the full optical train is calibrated, from the primary mirror to the de- 
tectors. 

• Unlike MEF, can be employed in the absence of sunspots (e.g., looking at plage regions 
or any other suitable source of polarization vectors). 

• Unlike MEF, can be applied to non-solar telescopes. 
There are also some drawbacks that should be mentioned: 

• Unlike MEF, requires design modifications. 

• Unlike BEM, requires an astronomical source of polarization. 

6. Conclusions 

The sub-aperture method proposed in this work provides a suitable means of calibrating a 
large telescope. This paper explores its possibilities and limitations by means of hound-and- 
hares tests. Potential advantages and disadvantages with respect to two other alternatives 
(the BEM and MEF) are discussed. 

The simulations indicate that the sub-aperture method is robust enough to meet the strin- 
gent ATST polarimetric requirements in practical situations. The amount of time required 
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Fig. 2. Elements of the error matrix E as a function of do from numerous mon- 
tecarlo hound-and-hares tests. The various lines represent different sources of 
polarization. Solid: Sunspot (slit position 1 in Fig 1). Dashed: Plage (slit posi- 
tion 2 in Fig 1). Dash-dots: Network element in the quiet Sun (slit position 3 in 
Fig 1). The horizontal dotted line represents the ATST polarimetric accuracy 
requirements. 
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Fig. 3. Cross-talk from I to V as a function of the sub-aperture calibration 
accuracy (r] t = 0). Solid: No pointing error (r] x ). Dashed: 77^=0.75 arc-seconds. 
Dash-dots: rj x =1.5 arc-seconds. The horizontal dotted line represents the ATST 
polarimetric accuracy requirement. 
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for the calibration procedure is probably inferior to ~20 minutes. Most of the time is spent 
observing the various configurations of the calibration optics, which is the same procedure 
employed at the DST (Sacramento Peak) or the VTT (Tenerife) for the polarimeter cali- 
bration. This is done currently on a routine basis about once per day during the observing 
runs. 
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